1 



Resonances in three-body systems with short and long-range interactions 

O ■ E. Garrido a , D.V. Fedorov b , and A.S. Jensen b 
O 



S3 . 

b 



in 



> 
o 



Instituto de Estructura de la Materia, CSIC, Serrano 123, E-28006 Madrid, Spain 



Department of Physics and Astronomy, University of Aarhus, DK-8000 Aarhus C, 
Denmark 



The complex scaling method permits calculations of few-body resonances with the 
correct asymptotic behaviour using a simple box boundary condition at a sufficiently large 
distance. This is also valid for systems involving more than one charged particle. We first 
apply the method on two-body systems. Three-body systems are then investigated by 
q ; use of the (complex scaled) hyperspheric adiabatic expansion method. The case of the 
2 + resonance in 6 Be and 6 Li is considered. Radial wave functions are obtained showing 
the correct asymptotic behaviour at intermediate values of the hyperradii, where wave 
functions can be computed fully numerically. 

3 i 1. INTRODUCTION 

c 

Use of the complex energy method permits to define resonances as poles of the S- 
matrix in the lower half of the momentum plane. The asymptotic behavior of the wave 
function is then known analytically as given by the Hankel function of first kind when only 
short-range interactions are involved, or by some combination of the regular and irregular 
Coulomb functions when the Coulomb potential enters. Resonance wave functions are 
then usually computed by imposing the correct analytic asymptotic behaviour. 

However the problem arises when dealing with more than two particles and when the 
long-range Coulomb interaction is involved. In this case the effective charge is not well 
defined, and therefore the correct asymptotics to be imposed to the solution is not known. 

The purpose of this work is to show that even in this case accurate resonance wave func- 
tions with the correct asymptotics can be computed. This is done by use of the complex 
scaling method, which transforms the exponentially divergent resonance wave function 
into a function that falls off exponentially, exactly as bound states do. This fact permits 
computations of resonances using the same numerical techniques as for bound states. In 
particular, a simple box boundary condition should be enough to obtain resonances after 
complex scaling. Nothing is then needed to be known and imposed in advance about the 
asymptotics, but still the correct asymptotic behaviour is recovered. 

In the following section we briefly describe the basis of the method. In section 3 we test 
it in the two-body case, for which the correct asymptotics is known also for two charged 
particles, and it can then be directly compared to the numerical solution. In section 4 
we apply the method for a three-body system, in particular for the 2 + resonances in 6 Be 
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(a+p+p) and 6 Li (a+p+n). The convergence and asymptotic behaviour of the solutions 
are discussed. We finish in section 5 with a short summary and the conclusions. 

2. THEORETICAL FORMULATION 

A system made of n constituents (n > 1) is usually described by a radial coordinate 
(relative distance for n—2, hyperradius for n—3, ■ ■ •) and 3n — 4 angles. Use of the complex 
energy method permits to obtain the resonances of such a system as poles of the S'-matrix 
in the lower half of the momentum plane. This implies that the radial wave functions of 
the resonances behave, for short-range potentials, as: 

/ reg ( r ) ^ ^H { l\nr) ™ e \K\rzine Re i(\K\rco 8 e R -<;n/2-*/4) (1) 

where is the Hankel function of first kind and order £, and k = ^2mE/h, where m 
is some normalization mass and E = \E\e~ l2eR is the n-body resonance energy. When 
the Coulomb interaction is involved, the Hankel function is replaced by the regular and 
irregular Coulomb functions, such that: 

/rcs(r) -> Fzfa «r) - iG^T], Kr) ™ e Mr sin 9 Re i(\n\r cos 9 R ~r,ln(2 K r)) ^ 

The first problem to face when computing resonance wave functions is the one coming 
from the exponential divergence observed in Eqs.flT]) and ([2]). For complex energy, ordi- 
nary continuum wave functions also diverge, and the distinction between them and the 
resonance wave functions is not trivial. 

This problem is solved by the complex scaling method, where the radial coordinates 
are rotated into the complex plane by an arbitrary angle 9 (r — > re 10 ). After this trans- 
formation, the radial resonance wave functions, except for some oscillatory term, behave 
asymptotically like f Ies {r) — ^ e _ ' K ' rsin ^ _e ' fl \ both for short and long-range potentials. 
Thus, as soon as 9 > 9r the resonance wave functions behave as bound states, and they 
can be computed following the same numerical procedures. In particular, since the radial 
wave functions go exponentially to zero, it is possible to obtain them with a simple box 
boundary condition f(r max ) = 0, with the only requirement that r max is large enough such 
that the wave function at that distance is orders of magnitude smaller than its maximum 
value. A box boundary condition is known to work for bound states, and should also 
work for resonances after a complex scaling transformation. 

An additional problem arises because for more than two particles the Coulomb charge 
T) (and to some extent also the index £) must be computed numerically. 



3. TWO-BODY CASE 

A two-body system is a good test of the method, since for this particular case the 
correct asymptotics is known also when the Coulomb interaction is present. 

Assuming a central two-body potential V(r), after complex scaling, the differential 
equation to be solved is simply: 

|L + | Up _ e *V (re ») _ /W(r) = (3) 



3 




Figure 1. Thick curves: Real and imaginary parts of the complex rotated radial wave 
function of the 2~ resonance in 16 F obtained with a box boundary condition. Thin curves: 
Correct asymptotics of the radial wave function as given by Eq.([2]). 



where /i is the reduced mass, and i is the relative orbital angular momentum. 

As an example we take a two-body system made by an 15 core and a proton. We 
use a nuclear interaction as given by the gaussian <i-wave potential in table 1 of [ [T] 
plus, of course, the Coulomb repulsion. We then solve the previous equation by imposing 
f(r max )—0, with r maa .=40 fm. A solution is found at the complex energy E = 0.5 — z0.004 
MeV, that matches well with the experimental 2~ resonance in 16 F [ [2] . 

In Fig{T]the thick-solid and thick-dashed curves show the real and imaginary parts of 
the complex rotated radial wave function of the resonance. As mentioned above, in this 
case the asymptotics is well known, as is given by Eq.© with f = t and r] = fZ ± Z 2 a, 

where Z 1 and Z 2 are the charges, a is the fine structure constant, and p = ^2p,E/h. This 
asymptotic behaviour is shown in FigJT] by the corresponding thin curves. We can see 
that the numerical radial wave functions, obtained with a simple box boundary condition, 
reproduce the correct asymptotic behaviour. 

4. THREE-BODY CASE 

After testing the method for two-body systems, we now use the same procedure for 
three-body systems, for which the wave functions are computed by using the hyperspheric 
adiabatic expansion method [[3]. In this method the different radial wave functions are 
obtained by solving the coupled set of differential equations: 

("^ - ^ + XM p2 + ^ - Q nn ) f n (p)+ £ (-2P nn ,( P )^- p - Qnn'(p)) Up) = 0(4) 

where p is the hyperradius, m is a normalization mass, and P nn i and Q nn ' are the non- 
adiabatic terms coupling the different radial wave functions [ |3] . 

In the previous equations the key quantities are the effective potentials A n (p) that are 
obtained as the eigenvalues of the angular part of the Faddeev equations [ [3] . Accurate 
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calculation of these eigenvalues is the first and essential requirement needed to obtain 
reliable radial wave functions. 

4.1. 2 + resonances in 6 Be and 6 Li 

The main properties of 6 He are accurately described when treating this nucleus as a 
three-body system made by an a-particle and two neutrons. Since both the alpha-neutron 
and the neutron-neutron interactions are well known, this nucleus appears as an almost 
perfect test for all the available numerical three-body methods. 

For the same reason, 6 Be (a+p+p) and 6 Li (a+n+p) are specially appropriate to in- 
vestigate three-body systems when more than one particle is charged. In particular, it is 
well established that 6 He has a 2 + resonance with energy and width (£' R ,r B ) = (0.83,0.11) 
MeV. The same analog 2 + state has been found in 6 Li and 6 Be at (1.67,0.54) MeV and 
(3.04,1.16) MeV, respectively [ 3]. Realistic detailed calculations using the (complex 
scaled) adiabatic expansion method concerning the 2 + resonance in 6 He can be found in 

[EE]. 
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Figure 2. Converged A-functions (see Eq.(j4])) Figure 3. The same as in Figf5]for the 2 + 
included in the calculation of the 2 + reso- resonance in 6 Li. 
nance in 6 Be. 

In Figs. [2] and [3] we show the A-functions used in the calculation of the 2 + resonance 
in 6 Be and 6 Li, respectively. Partial waves with £ x and l y up to 10 have been included. 
A maximum value of the hypermomentum [K — 2n + £ x + ly) equal to 20 has been used 
for all the components, but for the most relevant ones the maximum value of K has been 
increased up to 200 for s- waves, 90 for p-waves, and 60 for (i-waves. In total, slightly more 
than 2500 hyperspherical harmonics have been included in the calculation. The computed 
A's have then converged at least up to 100 fm. A reduction by a factor of 2 in the basis 
size produces A-functions indistinguishable from the ones shown in the figures. 

It is important to note that 6 Li, contrary to 6 Be and 6 He, is not borromean (the neutron 
and the proton can bind into deuteron). This is actually revealed by the lowest A in Figj3l 
that diverges parabolically to — oo. Therefore, the calculation for 6 Li requires additional 
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components that are forbidden in 6 Be and 6 He. These components correspond to the ones 
with zero isospin in the neutron-proton channel. 

Once the A-functions have been computed, the remaining step is to solve the coupled 
set of differential equations (j3|. According to the discussion in the previous section, 
after complex scaling a box boundary condition f n {pmax)=Q should be enough to obtain 
accurate three-body resonance wave functions. However, as shown in Figsj2] and [3] the 
effective potentials entering in Eq.pj have been accurately computed only up to p max =100 
fm, that is too little to expect the box boundary condition to work. Furthermore, accurate 
computation of the A's up to a p max value several times larger is too expensive from the 
computation time point of view. 
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Figure 4. Real parts of the radial resonance Figure 5. The same as in FigJUfor the 2" 
wave functions of the 2 + resonance in 6 Be. resonance in 6 Li. 



It is not difficult to see [ [3] that at large distances the A-functions go linearly with 
p (this is actually obvious from Figs. [2] and [3]). Also the Q nn functions entering in 
Eq.fllJ go like 1/p 2 , while P nn /,and Q nn i {n ^ n') go to zero faster. We have then, for 
p larger than 100 fm, used extrapolations of the A's, P's and Q's in Eq.(jl]) according to 
A n (p) = A Xn p + B Xn + C Xn /p, Q nn {p) = A Qn Jp 2 + B Qn Jp 3 + C Q Jp\ and as A PiQ /p 3 + 
Bp,q/ p 4 + Cp,q/p 5 for the non-diagonal P's and Q's. We have then solved the coupled 
equations (jl]) using the numerical A's in Figs. [2]and[3]up to 100 fm, and the extrapolations 
for larger p's. In Figs. H] and [5] we show the real parts of the complex scaled radial wave 
functions obtained in this way for the 2 + resonance in 6 Be and 6 Li, respectively. For 6 Be 
we use a complex scaling angle of 0.15 rads, and p ma s=700 fm is enough to obtain a 2 + 
resonance with energy and width (Er,Tji) = (2.94,1A5) MeV, that agree rather well with 
the experimental values [H]. For 6 Li the complex scaling angle is 0.10 rads and p max =1000 
fm. We obtain (P^,r^) = (1.67,0.50) MeV, that also agrees well with the experiment [H]. 
From the figures it is now clear that 100 fm is certainly not large enough to impose a box 
boundary condition at that distance. 

In the extrapolated region (p > 100 fm) the effective potentials entering in Eq.(j3]) go 
like A/p + B/p 2 + ■ • ■, where the A and B coefficients are computed numerically. For 
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Figure 6. First computed radial wave func- 
tion (thick curves) for the 2 + resonance in 
6 Be and the corresponding asymptotics (thin 
curves) . 
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Figure 7. The same as in FigJH] for the 2 + 
resonance in 6 Li. 



such a kind of potential the solutions are known to go asymptotically as in Eq.([2]), where 
£ and i] are easily related to A and B. In Figs. [6] and [7| we show the first radial wave 
function for the 2 + resonance in 6 Be and 6 Li, respectively. The thick curves are the 
numerical solutions, while the thin curves show the analytic asymptotic behaviour as 
given by Eq.([2]). As seen in both figures the matching with the numerical calculations 
is very good in the "extrapolated" region, as expected. The remarkable fact is that this 
agreement is also excellent already around 50 fm, in the region where the purely numerical 
effective potentials are used. 

In Figs. [8] and [9] we observe how the same behaviour appears for the second radial 
wave function. Similar results are also found for the other radial functions. Therefore the 
numerical resonance radial wave functions obtained with the converged effective potentials 
shown in Figs. [2]and[3], are consistent with the extrapolations used for the A's, P's and Q J s 
in Eq.(jl]), and also these solutions have converged to the expected asymptotic behaviour 
already at p- values around 50 fm, in the region where full accurate numerical calculations 
can be performed. 

5. SUMMARY AND CONCLUSIONS 

After a complex scaling transformation resonance wave functions behave asymptotically 
as bound states, i.e. falls off exponentially. In this work we exploit this fact to obtain 
resonances simply by using a box boundary condition at a sufficiently large distance. In 
this way prior knowledge of the correct asymptotic behaviour of the wave function is not 
required, permitting then to compute resonances for systems for which this asymptotics is 
not known, in particular for three-body systems involving more than one charged particle. 

After testing the two-body case (for which the correct asymptotics is known), we have 
investigated the case of the 2 + resonance in 6 Be and 6 Li. These two nuclei are specially 
appropriate, since the uncertainties coming from the two-body interactions are small. 
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Figure 8. Second computed radial wave func- Figure 9. The same as in FigJHlfor the 2" 
tion (thick curves) for the 2 + resonance in resonance in 6 Li. 
6 Be and the corresponding asymptotics (thin 
curves) . 



The (complex scaled) hyperspheric adiabatic expansion method is used. Accurate and 
converged effective potentials are obtained up to p ma x= 100 fm, and beyond this distance 
extrapolations are used. These extrapolations permit to know the correct asymptotic 
behaviour, that should be the right one at least for distances much larger than p max . 

When solving the radial equations with a box boundary condition at a large value of p 
(700 fm for 6 Be and 1000 fm for 6 Li) we have found that at rather modest distances (p ~50 
fm) the radial wave functions already match with the asymptotic behaviour obtained from 
the extrapolated effective potentials. This means that the numerical effective potentials 
obtained for p < p max are consistent with the extrapolations used to solve the radial 
equations. Furthermore, use of these numerical potentials is enough to obtain radial 
resonance wave functions that have already reached the asymptotic behaviour. This 
implies that observables related to the asymptotics of the wave functions can be safely 
computed using the resonance wave functions obtained by this procedure. An example is 
the energy distributions of the fragments after decay of the resonance. The experimental 
energy distributions are related to the energy distributions of the fragments at distances 
where the correct asymptotics has been reached [[7j. 
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